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Motivation 


• Success of output error estimation and 
adaptive mesh refinement in goal- 
oriented simulations 

- Automatic and user-independent 
production databases 

• Challenges of simulation-based design 

- High CFD expertise 

► Reliable mesh generation, long setup time 

► High cost due to repeated evaluation of 
objectives on fine, hand-crafted meshes 
or high uncertainty due to inappropriate 
meshes 



Objectives 


Adaptive discretization 
of aerodynamic shape optimization problems 


Accuracy 


a r 


Automation 


• Improve design confidence 


• Reduce level of CFD expertise 
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- Direct control over objective 
function discretization error 


- Eliminate the requirement to hand-craft 
general meshes appropriate for all 
candidate designs 
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- Shorten problem setup time 
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Progress toward improved efficiency 

• Reduce cost by systematically increasing the depth of refinement as the 
design improves 

- Progressive optimization strategy 

- Investigate challenges of dynamic error control 


Previous Work - Infrastructure 


1 . Embedded-boundary Cartesian mesh method 

• Arbitrarily complex domains, efficient and accurate 

• Irregularity confined to body intersecting cells 

2. Incremental strategy for h-refinement 

Adjoints 



3. Aerodynamic shape optimization 

• Gradient computation 









See AIAA Paper 
2013-0543 
(Smith et al.) for 
applications 


4. Output error estimates 

• Adaptive mesh refinement 



Optimize 
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Optimization: Discrete Formulation 


Modify Geometry 


Mesh & Flow Solve 


Evaluate Objective 


Compute Gradient 


Optimize 


min J (X, Q) 

A 

subject to 

R(X,Q) = 0 VIefi 

• Steady Euler equations 

• Gradient-based optimization 

- BFGS 

- SNOPT 

• Shape optimization 

M = f[T(X)] 



d J 
dX 
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Role of Adjoints 


Gradients 

J = f(x, Q) 

e-g- Cd + (Cl 

dJ _ dJ d.J dQ 
dX ~ dX + <9Q dX 

, 9R 9R dQ 

0 = 1 — 

dX dQ dX 


H 


dJ _ dJ 
d X ~ dX 







Role of Adioints 




Moo — 1 .1 , oc — “25° 

J = Cn + 0.2 Ca 


Density Adjoint 


Addition of mass 
decreases functional 


Not sensitive 


Addition of mass 
increases functional 


• Control problem 


-Optimal shape design: adjust design variables to control the flow and 
improve performance 


-Error analysis: adjust mesh refinement to control discretization errors 
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Linearization Details 


• Objective function gradient 


d J 
dX 


dj 

dX 
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dM 

~dT 



• Mesh sensitivities: infinitesimal perturbations are confined to cutcells 



Triangle to cut-cell connectivity established on-the-fly as the design 
evolves: triangulation connectivity and topology allowed to change 


Error-Estimate Details 



J( Qh) « J(Qh) - (K) r R(Q") - Wfc - ^) X R(Q") 
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Bound on remaining error in each 
coarse cell k 
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ek = X] (^Q ^ l ) T r (Ql) 
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Net functional error E? = 

fc =0 


Remaining Error 
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Given a user specified tolerance TOL, refine until E < TOL 


In practice, specify number of cycles, mesh-growth factor per cycle and 
cell-budget 
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Basic Framework Integration 


Modify Geometry 


Initial Mesh 



Adapt & Solve 


Evaluate Objective 


Compute Gradient 


Optimize 
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• Integration into existing, fixed mesh, 
optimization framework 

- Build sequence of adapted meshes 

- Pass values of objective and gradient 
from finest mesh to optimizer 

Multilevel parallelism 

- Mesh sensitivities in 
stand-alone code 


» In each design iteration, perform fixed 
(user specified) number of adaptations 

- Fixed depth strategy 

- Very robust and precise control over 
computational resources 

- May be inefficient 



Basic Example 



• Demonstrate numerical optimization with adaptive meshing 


• Study mesh convergence of 
objective function, its error 
estimates and gradients 

• Find angle of attack to minimize 
drag coefficient 

- Transonic flow, IVL = 0.8 • 

- NACA 001 2 airfoil L J 

- J = Cd, X = a 


- Initial design: cx, = 2° 
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Mesh Setup 



• Demonstrate numerical optimization with adaptive meshing 


• Study mesh convergence of 
objective function, its error 
estimates and gradients 

• Fixed-depth strategy 

- 8 adaptive refinements at 
each design iteration 

- Initial mesh ~1 ,700 cells 

- Final mesh -25,000 cells 























































































Near- field view of initial mesh 
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Optimization Convergence History 
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-2 


-3 


-4 


-5 


-6 


• Optimizer minimizes drag in 7 iterations 

• Gradient reduced by almost 5 orders of magnitude 

• Angle of attack history: 2°, 1 °, -0.5°, 0.01 °, -0.001 0 


Gradient 


Final Meshes After 8 Adaptations 
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Output Mesh Convergence 




Mesh convergence of drag and gradient at selected design iterations 

• Drag and gradient are well converged on meshes with -10,000 cells 

• Sign predicted correctly even on the coarsest mesh 


Convergence of Error Estimates 




• Key parameter to safeguard oversolving and transfer optimization 
to next mesh 
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Mesh Efficiency of Fixed-Depth Strategy 



Search Direction 

Mesh too fine Mesh too coarse 

• Angle of attack history: 2°, 1 °, -0.5°, 0.01 °, -0.001 0 


Objectives in Quadratic Form 


• Frequently use objective functions that contain quadratic terms 

- Penalty terms, e.g. ( Cl - T ) 2 

- Inverse design, e.g. j = J (p - P target ) 2 dS 

• As working variable approaches its target, adjoint variables vanish 



• Consequences include vanishing error estimates as optimality is 
approached, which effectively terminate adaptation, as well as 
strongly non-monotone error convergence 


Quadratic Example 


• Find angle of attack to match a target lift coefficient 

- Transonic flow, Moo = 0.8 

- NACA 001 2 airfoil 

- J = (Ci - 0.55) 2 , X = ex 

- Initial design: cx,= 0° 

- Final design: cut « 2° 

• Fixed -depth strategy 

- 9 adaptive refinements at each design iteration 

- Initial mesh -1,700 cells; final mesh -35,000 cells 


Convergence Histories 





• Optimizer matches lift in 6 iterations 

• Error convergence satisfactory in early design iterations, but 
becomes non-monotone and errors vanish at optimality 
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“Companion” Functional 



• Use a companion functional to eliminate numerical artifacts 
for quadratic objectives 

- Objective function working variable is used for error 
control and drives adaptation 

- Objective function drives design 

• Possible to implement at no additional cost 

- Arrange computations to use error estimates from the 
penultimate adaptation cycle and solve objective 
function adjoint only on the finest mesh 
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Quadratic Example 



• Find angle of attack to match a target lift coefficient 

- J = (Ci- 0.55) 2 , X = a 

- Jec = Ci, Error Estimate = £ 

• Compute conservative error 
estimate in objective function 

J=({Ci±s)~ 0.55) 2 • 

J < (Ci- 0.55) 2 ± A 

A = |2 (Ci - 0.55)e| + £ 2 
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Error Estimate 


Quadratic Example 



• Eliminated numerical artifact of 
vanishing error estimate near 
optimality 


Error Estimate 



• Objective function error estimate 
smoothly decreasing in all design 
cycles 
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Sonic-Boom Mitigation 



Drive vehicle shape by prescribing quieter near-field signals 



Inverse Design Model Problem 

Problem Setup 



• Prescribe a target signature from a known shape and 
verify that the optimization can recover this solution 


• 1 0 design variables that control body radius 


• Moo = 1 .5 and cx = 0° 


Initial Shape 
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Inverse Design Model Problem 


Problem Setup 


• Inverse design formulation: J — / (jp — Pt^et) 2 dS ath/L = 2 

Po oJ 


• Error control functional: J E c = - 17 - / (p — Poo) 2 dS 

P00J 



• Consider two cases 

1 . Fixed-depth strategy: 7 adaptation cycles in each design iteration 

2. Progressive optimization 
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Inverse Design Model Problem 

Fixed -Depth Strategy 
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Inverse Design Model Problem 

Fixed -Depth Strategy 



7 Adaptations, ~650k cells 
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Objective Function, 


Progressive Optimization 



• Minimize number of design iterations performed on finest mesh 

• Allow the designs to advance as far as possible on each level 
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Progressive Optimization 



4 Adaptations, ~1 30k cells 



Initial, 4 Adaptations Final, 4 Adaptations 


0.03 



°‘° 3 0 50 100 150 200 250 300 

Distance along sensor 


Terminated due to design variable bound violation near nose 
Peak-to-peak signal reduced by over a factor of five, smooth aft body 
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Progressive Optimization 



5 Adaptations, ~230k cells 



Initial, 5 Adaptations 



Final, 5 Adaptations 



• Terminated due to design variable bound violation near nose 

• Smoother nose shape, finer scales not resolvable on the previous 
mesh 
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Progressive Optimization 



6 Adaptations, ~350k cells 



Initial, 6 Adaptations Final, 6 Adaptations 



• Most work performed on this mesh: cost is roughly half of fixed-depth 
example per design iteration 

• Target matched to plotting accuracy but tip shape different from target 
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Sonic-Boom Inverse Design 



7 Adaptations, ~650k cells 



Initial, 7 Adaptations Final, 7 Adaptations 



Matched target shape in 1 2 design iterations 

Roughly a factor of two faster than fixed-depth strategy 

Mesh largely unchanged, could we re-use the same mesh? 
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Summary and Future Work 

• Developed framework for gradient-based optimization with 
capability to perform adaptive meshing in each design iteration 

- Promising approach to enhance accuracy, efficiency and 
automation of simulation-based design 

• Preliminary investigation of dynamic error control 

- Eliminated numerical artifacts in error estimates for 
objective functions in quadratic form 

• Future work 

- Use of error estimates to limit oversolving 

- Transfer of Hessian matrix as the design moves from 
mesh to mesh 

- Mesh re-use from nearby designs 
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Questions 


marian.nemec@nasa.aov 

michael.aftosmis@nasa.gov 


